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Coating the carbon nanotubes: Geometry of incommensurate long-range ordered 

physisorbed monolayers 
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The structures of long-range ordered physisorbed monolayer on a carbon nanotube are examined. 
Geometrical and energetical constraints determining the order of such monolayers are discussed. A 
number of symmetrically different, strongly bound adsorbate structures is found for Xe adsorbates, 
some of which differ very little in energy. The presented results suggest that the atomically uniform 
coating of carbon nanotubes is possible and offer a clear visualization of such coatings. 
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• Carbon nanotubes have already found a number of interesting applications 1]. However, to truly assess their 
potential for various applications, a detailed description of interaction of carbon nanotubes with atoms and molecules 
of their surroundings is required. In this context, it has been suggested that coatingof carbon nanotubes with 
atoms may improve their properties that are important for various uses (see e.g. Refs. y and0). The experimental 
demonstration of coating of carbon nanotubes with crystalline sheets of V2O5 layer-like structures has been presented 
111 Ref. U However, it is not clear yet how to create fine, ordered coatings of monoatomic thickness. The long- 
wavelength vibrational properties of such coatings have been examined in the literature (see e.g. Ref. 0), but the 
coatings are not easily visualized since the confinement of atoms/molecules on a nanotube (cylindrical) surface imposes 
a specific ordering of such structures, as will be shown here. 

In this letter, we are interested in geometry and binding energies of long-range ordered monoatomically thick 
coatings of carbon nanotubes. We consider the capped nanotubes and thus, the monolayer coating is on the outer 
side of the nanotube (exohedral adsorption). The microscopic arrangement of the adsorbed layer of atoms depends on 
the details of both the adsorbate-adsorbate and adsorbate-substrate interactions. These interactions are not precisely 
known yet with the exception of physisorbed gases. A large amount of experimental and theoretical data exists for 
physisorption of gases on graphite surfaces Q , which enables one to use this information and reconstruct the relevant 
potentials that can then be applied to the present case of interest 0, HIS El- Even with the information on potentials, 
it is still not trivial to find the ordered adsorbate structures that minimize the free energy of the system. This is mainly 
(y-^ ' because there are many different configurations with very similar energy The adsorbate-substrate potential exhibits 
corrugation and in the case of very large corrugation amplitudes, the adsorbate overlayer will follow the periodicity of 
the nanotube, since the energy it looses due to the unfavorable positions of its neighbors is less than the one gained 
due to the fact that all the atoms sit in the minima of the substrate potential, i.e. the overlayer structure will be 
commensurate with the substrate. In the other extreme case in which the corrugation of the substrate is negligible, the 
adsorbate layer will order in such a way to maximize the binding energy resulting from the interadsorbate interactions, 
and the resulting structure will be incommensurate 0. The corrugation of graphite is in between the two extremes, 
and a number of different phases of gases physisorbed on the graphite surface (commensurate and incommensurate) 
has been reported in the literature (5J, depending also on the substrate temperature. In what follows, I shall assume 
that the adsorbate structure is indeed incommensurate with the corrugation of the carbon nanotube. In that case, 
the adsorbate-nanotube potential can be treated as smooth, i.e. one can neglect its corrugation when calculating the 
O ■ total energy of the structure. 

Assume for the moment that all of the adsorbate atoms are at the same distance from the nanotube axis, i.e. that 
their centers lie on a cylinder of radius R, R = R t + h, where R t is the nanotube radius and h is the height of 
the coating above the surface of the nanotube. It can be a priori expected (I shall not do so) that the adsorbate 
structure with minimum energy will locally have a nearly hexagonal order. It is nevertheless not completely obvious 
how to reconcile this local arrangement of atoms with the fact that all adsorbate sites on a cylinder surface must be 
equivalent. 

We define two vectors (ai and a.2) on a cylinder surface connecting nearest neighboring adsorbate atoms along two 
inequivalent directions. Any of the two vectors can be specified by a polar angle spanned by the vector's starting and 
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ending points, and by the projection of the vector on the cylinder axis (z axis), 

ai = (a,aa), 

a 2 = 09,6/3). (1) 

The vectors are defined on a cylinder surface i.e. each of them specifies one of the cylinder geodesies which are simple 
coil curves (see Fig. ^| . The lengths of the vectors are not necessarily the same and it is in principle possible that the 
nearest neighbor distances are not the same along the two geodesies. 

In order to obtain a long range ordered adsorbate structure, all the adsorbate position vectors r; should be linear 
combinations of the elementary vectors ai and aa, i.e. 

Ti = r it j = iai + ja 2 

= (Rcos(ia + j/3), R sin(ia + j/3), iaa + jbf3), (2) 

where i and j are integers, and in the last line of the above equation ti j are given in Cartesian coordinates. There is 
an additional condition that must be satisfied by the parameters of the structure. Starting from one crossing point of 
two geodesies, the next crossing point should be reached by integer numbers (to, n) of primitive vector translations 
along both geodesies. This yields 

ma — n/3 + 2n 7 

maa = nbft. (3) 

The adsorbate structure is thus defined by a, (3, a, and b which satisfy Eq. (J3J. 

The total energy of the adsorbate structure can be written as a sum of all adsorbate-adsorbate interactions ,v, and 
interaction of the adsorbates with the carbon nanotube, V (this neglects the energy of adsorbate zero-point motion). 
The energy per adsorbate atom is thus 

oo n— 1 

E b = J2' v ^- r ^) + V(r .o), (4) 

i— — oo j—O 

where r< j is given by Eq. J2J, and the prime on the sum indicates that the term with i = j = is excluded. 

We now look for the structure with minimum adsorbate energy, i.e. we minimize Eq. (JIJ with respect to its 
parameters. Explicitly, the minimization procedure consists in choosing the height of the overlayer h, a and a 
parameters of the first geodesic, integer numbers to and n that determine and b of the second geodesic [Eq. Q] , and 
calculating the energy of the resulting structure from Eq. Since we minimize the internal energy of the overlayer 
and not its free energy, our results are applicable at zero temperature only. 

The indexing of the structures with (to, n) is not unique because the integers to and n depend on the two nearest- 
neighbor geodesies one chooses (e.g. instead of two geodesies defined by ai and a 2 one could have chosen the two with 
— ai and a 2 , or ai and ai ± a 2 , presuming that ai ± a 2 still points to the nearest neighbor, see Fig. However, the 
representation of adsorbate structures of interest in terms of planar adsorbate arrangements folded onto a cylinder 
surface is possible and allows for a unique indexing. All strongly bound structures of interest exhibit quasi-hexagonal 
order of adsorbates, with all six nearest neighbor atoms nearly equally separated from the atom in the center of 
a hexagon. Nevertheless, these quasi-hexagons cannot in general be obtained by mathematical wrapping of the 
perfect hexagonal planar lattice of minimal energy. This is a consequence of the fact that the sum of inter-adsorbate 
interactions is different when the adsorbates are positioned on the surface of a cylinder (Cartesian distances change) , 
and that the whole structure must "fit" onto the cylinder. To specify the geometry of adsorbate coating it is sufficient 
to characterize two points in the planar, nearly hexagonal structure obtained by unwrapping the coating (the cutting 
preceding the unwrapping is parallel to the z-axis). The vector w connecting these two points can be thought of as 
being wrapped around the z = cross-section of the cylinder. The whole planar structure is assumed to follow this 
"wrapping" procedure. This is quite similar to the usual way of indexing of carbon nanotubes Q. The wrapping 
procedure is sketched in Fig. ^ 

The vector w can also be expressed as 

w = riiaj — na 2 , (5) 

where a' x and a' 2 are basis vectors of the planar, nearly hexagonal adsorbate structure which, when wrapped around 
the cylinder, produce the vectors ai and a 2 in Eq. JIJ, and the integers to and n correspond to those in Eq. ©. If 
one chooses a[ and a 2 in the way sketched in Fig. ^ (angle between a^ and a 2 greater than 7r/2), the indexing of 
coatings in terms of to and n [or wrapping vector in Eq. JSJ)] becomes unique and that is what we shall use to specify 
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(4,7) 


3.710 


-231.54 





(4,12) 


3.730 


-242.23 


1 


(1,9) 


3.705 


-231.41 


1 


(1,14) 


3.740 


-242.10 


2 


(5,6) 


3.700 


-231.29 


2 


(5,11) 


3.720 


-241.88 


3 


(3,8) 


3.740 


-231.06 


3 


(3,13) 


3.720 


-241.72 



TABLE I: Summary of the key characteristics of four lowest energy coatings with different symmetries for (5,5) and (10,10) 
single- walled carbon nanotubes. Wrapping vector (w), height above the carbon nanotube surface (h) and binding energy (Eb) 
are specified. 

the symmetries of the coatings. Note that the minus sign in Eq. (jSJ is a convention: one can also choose the plus sign 
and redefine vectors and a! 2l so that the angle between them is smaller than n/2. Alternatively, m and n numbers 
are those pertaining to ai and Si2 vectors chosen in such a way that the angle they span is maximum, while aa(3b > 
[Eq.©]. 

We now investigate the symmetries and binding energies of the adsorbate coatings as a function of the carbon 
nanotube radius. Adsorption of xenon on nanotube bundles has been experimentally studied (see Refs. mm 

and 

molecular dynamics simulations of Xe adsorption on carbon nanotube have led the authors of Ref. to conclude 
that the corrugation of the Xe-nanotube potential is not large enough to impose commensurate Xe structure. That 
is why the Xe monolayer coatings are considered in the following, although the approach presented here can be easily 
extended to other adsorbates which physisorb on nanotubes. 

The Xe-Xe (v) and Xe-nanotube (V) potentials are constructed as described in Ref. (see also Refs. I^ ITi)) . These 
are not the most precise potentials, but are fairly simple and shall suffice to illustrate the most important points of 
this article. 

Table 1 summarizes the results of calculations of binding energies per adsorbate atom for several lowest energy 
configurations of Xe coatings and for two different radii of a single wall carbon nanotube. The configurations are 
numbered as 0,1,2,... and are ordered such that their energy increases with the increase of the configuration number. 
The radii chosen correspond to (5,5) and (10,10) carbon nanotubes. The symmetries of the coatings are denoted. 
The four lowest energy Xe structures on (5,5) carbon nanotube are plotted in Fig. [5] The parameters for the lowest 
energy structure (m = 4, n = 7) on (5,5) nanotube are (a = 0.506 rad, (3 = —0.609 rad, a = —5.650 A/rad, b — 2.681 
A/rad), while for (4,12) coating on (10,10) nanotube, (a = 0.302 rad, (3 = -0.423 rad, a = -10.87 A/rad, b = 2.583 
A/rad). 

One should bear in mind that it is possible to find a configuration of Xe coating with the energy in between those 
specified in Table 1. Such a structure can be realized by slightly straining the intcradsorbate distances in the minimum 
energy configuration of a particular symmetry in such a way that the symmetry of the structure remains the same. 
For example, the minimum energy structure of a particular symmetry can be turned into the structure of the same 
symmetry but smaller binding energy, by applying forces to the coating along the z direction to slightly squeeze (or 
stretch) it. The discrete energies quoted in Table 1 correspond to the lowest energy configurations of a coating with 
specified symmetry and it is not possible to obtain lower energy by applying any set of infinitesimal changes to the 
adsorbate position vectors. Obviously, a large number of different symmetry coatings differing only a little in the 
binding energy is obtained. For example, for (10,10) nanotube highly bound coatings arc also those with wrapping 
indices (2,13), (6,10), (0,14), (8,9), (7,9) ... But the coatings with different symmetries do not all find the minimum 
energy configuration at the same distance from the tube center. This is due to the fact that it may be advantageous to 
the overlayer to slightly move away from the height at which it interacts with the nanotube most attractively, in order 
to "repack" itself in such a way to increase a part of the binding energy which is due to interadsorbate interactions. 
All the energies are in fact doubly degenerate, since for each choice of (a, /3, a, 6), there is a structure of exactly the 
same energy, but topologically different, with parameters (a, /3, —a, —b). This could be termed as helicity degeneracy. 
Futhermore, any rotation of the coating around the z-axis does not change its energy. 

Note also that the adsorbate coatings do not in general show periodicity in z -direction. This is due to the fact that 
the unwrapped 2D hexagonal lattice is not perfect (see Fig. Only when ma — 2/-7T, where / is a rational number, 
and m and a those from Eq. @, the adsorbate structure is periodic in z-direction. 

For multiwalled carbon nanotubes, the conclusions reached at here will not change significantly. In that case only 
the adsorbate-substrate potential change, and the procedure described above can applies equally well. Nevertheless, 
due to small energy differences between different configurations of the coatings, it is possible that the minimum 
energy configurations may exhibit different symmetries for single-walled nanotube and multiwalled nanotube with the 
same (external) radii. The presented approach can be also applied to the coatings which form inside the nanotubes. 
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Essentially the only difference in this case is that h is negative. The packing of fullerene molecules in carbon nanotubes 
has been considered in Ref. ITj. Packing of hard spheres in cylinders has been studied in Ref. [0| The authors of this 
reference emphasize the chirality aspect of the structures which is also obtained in this work. The lattices of hard 
spheres formed on cylindrical surfaces have been studied in Ref. Il7l 

There is an additional degree of freedom available to the coating which was effectively "frozen" in the presented 
approach. Namely, one could imagine a situation where the adsorbate monolayer additionally relaxes in such a way 
that its height above the nanotube becomes positional dependent. This is similar to the phenomenon of the overlayer 
"buckling" encountered in some adsorbate/planar surface systems of for some clean surfaces (see e.g. Ref. 18). The 
simplest form of "buckled" coating would be the one with the atomic centers positioned on an ellipse in the x — y 
plane (instead of a circle) . 

An interesting question which comes to mind, especially in view of such small energy differences among various 
configurations, is whether the inclusion of zero point motion energy would change the results presented in a qualitative 
way. The different structures are indeed expected to have inequal zero-point motion energies, since the potentials 
experienced by adsorbates differ in these configurations, at least slightly. This is probably a very small effect whose 
evaluation may require quite sophisticated techniques. 

In conclusion, we have considered geometrical and energetical constraints and requirements associated with coating 
of the nanotube with a monolayer of physisorbed atoms. We have theoretically predicted a number of different 
symmetry coatings, some of which differ only a little in their binding energy. The results should be of importance 
for future experiments directed towards functionalization of carbon nanotubes by such fine, monoatomically thick 
coatings. 

Stimulating and usefull discussions with D. Jurman, Dr. B. Gumhalter and Dr. H. Buljan are acknowledged. 
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FIG. 1: Nearly perfect hexagonal lattice of points. The way to turn this lattice into an adsorbate coating is illustrated by 
the dashed lines and "scissors" (cutting along the dashed lines and wrapping the w vector around z = cross-section of the 
nanotube). The basis vectors and a' 2 used to specify any wrapping vector w are denoted. Several minimum energy structures 
obtained for (5,5) nanotubes are indicated. Note that although all these structures are indicated on the same hexagonal lattice, 
they in fact belong to slightly different nearly perfect hexagonal lattices as discussed in the text. The inset displays ai and &2 
vectors and their components. 



FIG. 2: Four lowest energy Xe structures on (5,5) carbon nanotube. From left to right: (4,7), (1,9), (5,6) and (3,8) coatings. 



